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We provide a scan of the parameter space for neutralino-hadron scattering in the next-to-minimal supersym- 
metric standard model using an updated value for the strange quark sigma commutator These results also take 
into account constraints from WMAP data on the relic density and new constraints from the Large Hadron Col- 
lider We find that the resultant spin-independent cross sections are smaller in magnitude than those found in 
recent results obtained within the constrained minimal supersymmetric standard model, yet still great enough to 
feasibly allow for detection in the case of bino-like neutralinos. 
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It is now well established that a vast portion of our universe 
(almost a quarter) is comprised of a weakly interacting quan- 
tity known as dark matter. There is no satisfactory candidate 
within the standard model and thus its categorisation remains 
one of the great conundrums of modern physics. Currently, 
theories in which dark matter consists of slow-moving par- 
ticles are most favoured, on the basis of extensive evidence 
from rotation curves, galaxy clusters, gravitational lensing 
and anisotropies in the cosmic microwave background radi- 
ation [ ] . In particular, Weakly Interacting Massive Particles 
(WIMPs) are the leading candidates. There are many detec- 
tion experiments underway around the world that seek to es- 
tablish the presence of such WIMPs, both directly (includ- 
ing CDMS, DAMA/LIBRA, DRIFT, EDELWEISS, LUX, PI- 
CASSO, SIMPLE, WArP, XENON, XMASS and ZEPLIN- 
III [2- ]) and indirectly (including AMANDA, ANTARES, 
Fermi-LAT, IceCube and PAMELA [ - ]). It is vital to thor- 
oughly explore the properties of dark matter, particularly its 
interactions with hadronic matter, within a variety of different 
theoretical models, since this will help to guide direct efforts 
to detect it. 

Supersymmetric models have been quite popular within the 
past few decades - particularly the minimal supersymmetric 
standard model (MSSM). In these scenarios, dark matter takes 
the form of the lightest supersymmetric particle (LSP), usu- 
ally a neutralino; it is favoured as a dark matter candidate for 
its stable and weakly-interacting nature [ ]. Various authors 
in the past decade or so have used the MSSM and its variants 
to make predictions about its interaction with baryonic mat- 
ter via the calculation of spin independent neutralino-hadron 
scattering cross sections. However, in the light of recent 
results from the LHC, the MSSM is becoming increasingly 
hard-pressed to serve as a complete description of physics be- 
yond the standard model, because of the need to fine-tune its 
parameter space. The next- to-minimal supersymmetric stan- 
dard model (NMSSM), on the other hand, is under less threat 
[10], and is currently a model of considerable interest. 

In this Letter, we calculate the spin-independent cross sec- 
tion for neutralino-nucleon scattering within the constrained 
NMSSM, taking into account the constraints from the initial 
running of the LHC, as well as lattice QCD determinations of 
the light quark sigma commutators. We incorporate the latest 



relic density constraints from WMAP results and find a dras- 
tic reduction of regions in the NMSSM parameter space for 
which neutralino dark matter is viable. We show that the spin- 
independent cross sections fall within two main regions; the 
first corresponding to bino-like neutralinos with cross sections 
on the edge of detection limits, the second to singlino-like 
neutralinos with cross-sections far too small to be detected in 
current experiments. 

To begin, we briefly outline the relevant features of the 
NMSSM; for a comprehensive review, we recommend [11]. 
The MSSM introduces two neutral Higgs doublets H u and 
Hd to the standard model. The Higgs superfields contribute a 
Higgs mass term to the superpotential of the MSSM, 



W; 



MSSM 



= W Y 



nH u Hd, 



(1) 



where Wy represents the Yukawa couplings for the SM 
fermions and H u , H^ are the Higgs chiral superfields. In or- 
der to avoid extreme fine-tuning, it is necessary that the /i term 
and the scale of SUSY breaking both lie at the electroweak 
scale. It is unknown why the two scales should fall so close 
to each other (and far below the GUT scale) when \i itself has 
little to do with SUSY breaking; this is generally considered 
to be a problem of naturalness. 

Historically, the NMSSM was formulated as a convenient 
way of dealing with this "/i problem." In the NMSSM, /i is 
replaced by a gauge singlet chiral superfield S. An effective /i 
can thus be dynamically generated upon SUSY breaking, ex- 
plaining the coincidence of scales. This results in an expanded 
superpotential in the NMSSM, 



Wnmssm = Wy + \SH u H d 



K q3 
3 S ' 



(2) 



The promotion of ji to a singlet field does have some con- 
sequences for the neutralino-hadron cross-section, since it re- 
sults in a greater number of ways in which these particles can 
interact. Two extra Higgs fields are generated, such that the 
Higgs sector of the NMSSM consists of three neutral CP-even 
Higgs, two CP-odd Higgs and two charged Higgs. However, 
only the three CP-even Higgs are of relevance when formulat- 
ing neutralino-hadron spin-independent cross sections. 



The lightest neutralino itself may be described as a mixing 
of five neutral fields rather than the four of the MSSM, 



X = Z x iB 



Z x iW 



Z^zHa -\- Z^aHo 



Z^bS. 



(3) 



Here, B is the bino (superpartner of the U(l ) gauge field), W 
is the wino (superpartner of the W gauge field), Hi and H2 
are higgsinos (superpartners of the Higgs fields) and S is the 
singlino (superpartner of the singlet field). The behaviour of 
the neutralino as it interacts with hadronic matter is strongly 
dependent on its exact composition. A predominantly bino- 
like neutralino (99% bino), for instance, will be shown to yield 
a high spin-independent cross-section, relative to a singlino- 
like neutralino, since singlinos do not couple to sfermions, 
quarks or gauge fields. 

Before proceeding to the results, it may be helpful to outline 
the method used to calculate the spin-independent neutralino- 
hadron cross section, a si- (We refrain from addressing the 
spin-dependent component of the neutralino-hadron cross sec- 
tion, since it is typically several orders of magnitude below 
experimental sensitivity [12].) 




FIG. 1 : A neutralino-nucleon collision via Higgs exchange. 

Fig. 1 shows one example of an interaction between a neu- 
tralino (denoted by x) an d a nucleon (neutralinos and nucle- 
ons may also interact via squark exchange). The cross section 
for this kind of interaction involves matrix elements of the 
form (A\qq\A). The following gives a succinct expression for 
the contribution to a si of each quark [13-15]: 

(Nx\(X3qXXqq\xN) = a 3q (N\qq\N) 

fr q 

= ^ q m N , 

m q 

where we have used rriNfrq = (N\m q qq\N) for the sigma 
terms (the bottom half of Fig. 1) and a 3q encapsulates the 
relevant physics in terms of the amplitudes of each contribut- 
ing neutralino-quark interaction (the top half of Fig. 1). Sum- 
ming over light and heavy quarks gives a generic expression 
for a si = 4 (m^/ir) f 2 , where 
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m N 
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q=u,d,s 
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(4) 



This may be simplified further by noting that tun = 
(NlO^N) for a system at rest, where the trace of the energy- 
momentum tensor is given by 

^ = E m ^+E m ew-^<^<^. ^ 

q Q 

Taking the system to be at zero-momentum, such that 9fi = 
#q, this equation becomes 

m N = (N\ J2 m q qq\N) + (N\ ^ m Q QQ\N) 

q Q 

ley 

-(N\^G^G^\N). (6) 

The extra term involving G^G^ can be conveniently elim- 
inated by appealing to the Novikov-Shifman-Vainshtein [16] 
relation, which tells us that - (7a s / Sir) (N\G^ U G^\N) may 

be written as (7/2) £ Q m Q (N\QQ\N), 

m N = J2( N \™\ N ) + I X>Q<W|QQW 
q Q 

= Y,™NfTq + f m Q (N\QQ\N), 
q 

andthusm Q (A^|QQ|AT) = ^m N (l - ^ f Tq \ Hence, the 
final expression for / is 



a 3Q 



f V^ a SqfTq 2 f V^ 

m N ~ 2^ m + 27 Jtq 2^ m 

q—u,d,s H Q=c,b,t ^ 



/tq=(i- E M). 



(7) 



q=u,d,s 



The as q terms for the NMSSM themselves are derived by 
computing the amplitudes of the contributing Feynman dia- 
grams. In addition, the fr q terms have been derived from 
numerous studies in lattice QCD [17-19], which are updates 
of early estimates [20, 21]. 

In order to compute the evolution of the coefficients as q 
from the GUT scale to the EW scale, we modified mi- 
crOMEGAs, a code for calculating general dark matter prop- 
erties under supersymmetric physics, developed by Belanger 
et al. [ ]. In order to reduce the parameter space, we chose 
a constrained version of the NMSSM in which the scalar and 
gaugino masses were taken to be universal at the GUT scale. 
Thus, the free parameters are the universal scalar mass mo, 
universal gaugino mass rai/2, singlino trilinear coupling A K9 
Higgs-singlino trilinear coupling A\ (later referred to simply 
as A in this paper), tan/3 (the ratio of vacuum expectation val- 
ues of the neutral Higgses) and A. In addition, the effective 
Higgs mass \i was taken to be positive. One of the constraints 
imposed by micrOMEGAs (namely, the computation of the 
muon anomalous moment [23]) was relaxed. 

Initially, the light quark sigma term o\ was taken to be 47 
MeV, to correspond with a lattice determination of 07 = 47 ±9 
MeV, and the strangeness sigma term cr s was taken to be 
50 MeV, in accordance with a s = 50 ± 8 MeV from [17], 
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FIG. 2: Regions in the space of universal spin-1/2 and spin-0 masses allowed by relic density constraints 
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FIG. 3: Cross sections for the spin-independent neutralino-neucleon cross section for the parameter sets illustrated in Fig. 2a 

and Fig. 2b 



which was obtained by averaging two different lattice results. 
However, recent findings tend to favour an even lower value 
[24, 25]. Thus, these scans were repeated with cr s = 22 ± 6 
MeV and o\ = 47 ± 9 MeV, and the plots provided used these 
values. In addition, findings from WMAP and other obser- 
vations have placed constraints on the relic density Q to lie 
between 0.1053 and 0.1193 at 95% confidence level [26, 27]. 
This constraint places tight restrictions on the allowed pa- 
rameter space, where for a fixed A, the allowed regions are 
reduced to thin strips or lines in the (mo,mi/ 2 ) plane. Fi- 
nally, recent data from the CMS collaboration [ ] was used 
to place a lower bound on (mo^^), with similar results 
having been obtained from ATLAS [29]. Although this bound 



was originally formulated within the context of the MSSM, 
the spectrum of superpartners is quite similar in the NMSSM 
within this region of parameter space, so the bound still rep- 
resents a very good approximation in the present case. 

Sweeps of mo and mi / 2 were carried out for various values 
of A at high tan/3 (50), with A and A K fixed at 0.01 and -40 
respectively. Points that are allowed by both LHC and relic 
density constraints are plotted in Figure 2a. Noticeably, as 
rai/2 is increased the character of the neutralino changes quite 
significantly. For lower m 1 / 2 along each strip, the content 
is predominantly bino, and cross-sections are of the order of 
10 -9 pb. However, as the points cross the bino-singlino line, 
as illustrated in Figure 2a, the neutralino is almost entirely 



singlino and the cross-section sharply drops by several orders 
of magnitude. For this reason, the region immediately to the 
right of the bino-singlino line is most favourable as an area 
of interest, since neutralinos in the singlino-like region have 
interaction cross sections that are far too small to allow for 
detection in current searches (see Figure 3a). 

In addition, lower values of tan/3 (5 and 10) also yield al- 
lowed regions, although the highest cross-sections are found 
for tan/3 = 50 (see Figure 3b). However, in the low tan/3 region 
we do not observe any singlino-like behaviour (as evidenced 
in Figure 2b), since a limit is quickly reached beyond which 
the LSP is no longer a neutralino. Very low values of tan/3 
proved to be unfavourable, because of the presence of a Lan- 
dau pole. 

TABLE I: Example breakdown of quark flavour contributions 
at high tan/3 with a s = 22 MeV, o x = 47 MeV. 



Model 


q 


OL3q/m q 


fllh 


tan/3 = 50, A = -575 


u 


-1.179 x 10~ y 


0.0196 


rao = 436, m 1 / 2 = 510 


d 


-1.090 x 10~ 8 


0.1820 


A = 0.01, A* = -40 


c 


-1.179 x 10~ 9 


0.0538 


a S i = 8.678 x lO" 10 pb 


s 


-1.090 x 10~ 8 


0.1700 


ai = 47 MeV 


t 


-1.174 x 10~ 9 


0.0536 


a s = 22 MeV 


b 


-1.142 x 10~ 8 


0.5213 



Table I shows that in the NMSSM at high tan/3, the cross- 
section is dominated by the down-type quarks (particularly the 
bottom and strange quarks). This is similar to the finding in 



the CMSSM, although in the present case the bottom quark 
dominates by an even greater percentage. This is shown also 
for Table I with an updated cr s value; lowering cr s from 47 
MeV to 22 MeV seems to have the effect of reducing the cross 
sections by typically 30%. 

From these results, it seems that the constrained NMSSM 
does indeed produce a number of viable dark matter candi- 
dates, a si for the neutralino-hadron collision in this scheme 
is very strongly dependent on the composition of the neu- 
tralino itself. In the region where the neutralino is predom- 
inantly bino-like, a si is comparable to the values found in the 
CMSSM and within the reach of direct detection experiments. 
On the other hand, there is a second region in which the neu- 
tralino is predominantly singlino-like, where a si is negligibly 
small. Furthermore, given the sharp drop in the a si for these 
singlino-like neutralinos, our results suggest a possible sce- 
nario in which a discovery at the LHC may be compatible with 
a null result in direct detection dark matter searches. Further 
investigation will be necessary for other variations of super- 
symmetric models as new data from the LHC is produced. 
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